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Abstract 

This  work  presents  a  study  of  the  determination  of  uncertainties  in 
[UTC  —  UTC{k)}  needed  for  publication  in  the  Bureau  Inernational  des 
Poids  et  Mesures’s  (BIPM’s)  Circular  T  and  the  Key  Comparison  Database, 
as  required  by  the  Mutual  Recognition  Arrangement.  In  the  hrst  part  of  the 
paper,  an  analytical  solution  based  on  the  law  of  the  propagation  of 
uncertainty  is  derived.  In  the  second  part,  the  solution  is  verihed  numerically 
using  the  software  used  by  the  BIPM  for  the  generation  of  UTC. 

(Some  hgures  in  this  article  are  in  colour  only  in  the  electronic  version) 


1.  Introduction 

Coordinated  Universal  Time  (UTC),  the  worldwide  time 
standard,  is  disseminated  monthly  through  publication  of 
[UTC  —  UTCik)]  in  Circular  T  by  the  Bureau  Inernational 
des  Poids  et  Mesures  (BIPM).  UTC(/:)  is  a  realization  of  UTC 
by  laboratory  k.  These  values  before  2005  have  been  published 
without  their  uncertainties,  but  the  rules  of  the  CIPM  Mutual 
Recognition  Arrangement  (CIPM  MRA)  require  the  evaluation 
of  this  uncertainty.  This  paper  reports  the  first  steps  towards 
computing  these  uncertainties. 

UTC  is  derived  from  International  Atomic  Time  (TAI)  by 
the  addition  of  leap  seconds;  while  TAI  is  derived  from  the 
Free  Atomic  Timescale  (EAL)  through  the  addition  of  pre¬ 
announced  frequency  steers  determined  by  comparison  with 
a  weighted  set  of  primary  frequency  standards.  EAL  is  a 
worldwide  weighted  average  of  a  large  number  of  free-running, 
effectively  uncalibrated,  frequency  standards  [1-3].  The 
uncertainty  in  the  determination  of  EAL,  TAI  and  UTC,  as  steps 
in  the  realization  of  Terrestrial  Time  (TT),  is  affected  by  three 
major  elements:  clock  variations,  the  means  of  comparisons 
of  remote  clocks  (time  transfer)  and  the  time-scale  algorithm. 
The  uncertainties  of  time  transfer  are  particularly  significant 
over  averaging  times  of  up  to  a  few  tens  of  days  and  also 
influence  the  uncertainty  in  the  access  of  the  participating 
laboratories  k  to  UTC  (in  other  words,  the  uncertainty  of 
[UTC  -  UTC{k)\). 


In  this  work,  a  study  of  the  determination  of  the 
uncertainties  of  [UTC  —  UTC{k)]  is  presented.  In  the  first 
part  of  the  paper,  an  analytical  solution  based  on  the  law  of 
the  propagation  of  uncertainty  is  derived.  The  second  part 
presents  a  numerical  verification  of  the  analytical  results  using 
the  software  used  for  the  generation  of  UTC. 

2.  The  algorithms  of  EAL,  TAI  and  UTC  and  the 
working  hypothesis 

The  uncertainty  of  [UTC  —  UTC{k)]  can  be  derived  using  the 
general  equation  of  the  Eree  Atomic  Timescale  (EAL),  which 
is  defined  using  the  ALGOS  algorithm  [1-3]  as 

N 

EAL{t)  =  Y,^i\.hii.t)  +  h\{t)],  (1) 

i  =  \ 

where  N  is  the  number  of  the  atomic  clocks,  wt  the  relative 
weight  of  the  clock  H/,  hi(t)  is  the  reading  of  clock  H/  at 
time  t  and  h'-(t)  is  the  prediction  of  the  reading  of  clock  H/ 
to  guarantee  the  continuity  of  the  time  scale.  The  weight 
attributed  to  a  given  clock  reflects  its  long-term  stability,  since 
the  objective  is  to  obtain  a  weighted  average  that  is  more  stable 
in  the  long  term  than  any  of  the  contributing  elements  [4, 5]. 
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The  evaluation  of  uncertainties  in  [UTC  —  UTC{k)] 


The  weights  of  the  clocks  obey  the  relation 

N 

=  (2) 

i=l 

Subtracting  the  same  quantity  from  both  sides  of  equation  (1), 

N  N  N 

EAL{t)  —  ^  Wihi(t)  =  ^  Wi[hi{t)  +  ^  Wihi(t). 

i  =  \  i  =  \  i  =  \ 

Using  equation  (2)  and  rearranging, 

N  N 

Y,WiiEAL{t)  -  hi{t))  =  J^wth',(t).  (3) 

i=l  i=l 

Setting 

Xi(t)  =  EAL(t)-hi(t),  (4) 

it  is  clear  that  equation  (3)  is  of  the  form 

N  N 

wiXiit)  =  Wih'iit).  (5) 

i  =  l  i  =  l 

The  software  package  termed  ALGOS  is  used  in  the  Time 
Section  of  the  BIPM  to  generate  UTC.  Weights  are  determined 
from  the  variance  of  monthly  average  frequencies,  subject  to  a 
maximum  value  [5].  The  data  used  by  ALGOS  take  the  form 
of  the  time  differences  between  readings  of  clocks,  written  as 

=  hj(t)  -  (6) 

Equation  (5)  in  conjunction  with  the  A  —  1  equations  (6)  results 
in  a  system  with  N  equations  and  N  unknowns: 

I mxi(t)  =  Wih'iit) 

\xi(t)  -  Xjit)  =  Xijit) 

The  solution  is 

N 

Xj(t)  =  EAL(t)  -  hj(t)  =  wdh'iit)  -  Xij(t)l  (8) 

i=l 

The  difference  between  any  clock  H ^  and  EAL  (8)  depends 
on  weights,  clock  prediction  and  measured  clock  differences. 
The  clock  Hj  may  also  represent  a  UTC(7)  time  scale; 
therefore,  xj  (t)  may  also  be  interpreted  as 

Xj  =  EAL  -  UTC(j) 

having  dropped  the  time  instant  t  in  the  notation  for  simplicity. 

The  predictions  and  the  weights  are  fixed  by  appropriate 
algorithms  based  on  past  clock  behaviour,  and,  in  equation  (8), 
they  can  be  considered  as  time-varying  deterministic 
parameters.  Suboptimal  estimation  of  these  parameters  would 
affect  the  uncertainty  of  TAI  as  a  realization  of  the  TT,  but  they 
do  not  affect  the  knowledge  of  the  difference  between  EAL  and 
clock  Kj ;  the  measures  xtj  are  thus  the  only  contributors  to  the 
uncertainties  in  xj .  The  contribution  of  the  uncertainty  given 
by  measures  of  clocks  located  inside  the  same  laboratory  is 
considered  negligible,  and  in  terms  of  the  ALGOS  algorithm 
it  can  be  ignored  because  it  is  inseparable  from  the  noise  of  the 
clock. 


Because  [UTC  —  EAL]  depends  only  on  pre-determined 
leap  seconds  and  frequency  steers  that  do  not  add  uncertainty, 
and  because  UTC(^)  is  deternfinistically  derived  from  the 
laboratory  clock  readings,  the  uncertainties  of  [UTC  — 
UTC(k)]  are  identical  to  the  uncertainties  of  [TAI  —  UTC(k)], 
[EAL  —  UTC(k)]  and  [EAL  —  H^].  Thus,  the  uncertainties 
of  the  links  between  laboratories  are  the  only  source  of  the 
uncertainty  of  [UTC  —  UTC(k)].  In  the  following  sections  its 
propagation  will  be  studied. 

2.7.  The  law  of  propagation  of  uncertainty 

According  to  [6],  the  uncertainty  in  the  Xj(t)  can  be  found 
using  the  law  of  the  propagation  of  uncertainty.  Defining  y 
as  a  generic  quantity  indirectly  measured  by  means  of  direct 
measurements  of  the  input  quantity  x/ , 


y  =  /(xi,X2, 


The  expression  of  the  law  of  the  propagation  of  uncertainty 
[6]  is  given  by 


u 


2 

V 


M-\  M 


i=l  k=i+\ 


dxi 


(9) 


where  the  first  term  corresponds  to  the  effect  of  the 
uncertainties  on  the  input  quantities  x/  and  the  second  term 
accounts  for  the  correlation  between  them. 

The  uncertainty  of  the  quantity  Xk  =  [EAL  —  UTC(k)] 
defined  in  equation  (8)  plays  here  the  role  of  the  indirect 
quantity  y  and  the  uncertainty  contributions  are  only  due  to  the 
measurement  noise  of  the  links  Xij(t).  Applying  equation  (9) 
to  the  model  (8)  yields 


.V  /  ^  X  z  N-l  N 

=  E(^)  <-2EE 

/  =  1  /  =  1  k=M 


dxj  dxj 
dxij  dxk, 


_  ^iXij,Xkj) 


=  E  2  E  E  (10) 

i=l  i=l  k=i+l 

where  «2.  = 

The  weights  of  the  clocks  are  available  from  the  BIPM 
website,  and  the  uncertainties  of  links  between  the  clocks  [7] 
are  published  in  Circular  T  (see  the  BIPM  Time  Section’s  ETP 
server  at  http  ://www.bipm. org/j  sp/en/TimePtp .j  sp) . 

The  propagation  of  uncertainty  (10)  could  as  well  be 
expressed  by  a  matrix  formulation  using  a  multivariate 
weighted  approach  [8,9];  the  scalar  approach  presented  in 
this  paper  allows  one  to  clearly  understand  how  the  different 
components  contribute  to  the  combined  uncertainty. 

It  can  be  demonstrated  that  the  ALGOS  algorithm  would 
generate  the  same  results  if  each  laboratory’s  clocks  were 
replaced  by  a  single  ‘equivalent’  clock  whose  reading  was  the 
weighted  average  of  the  individual  clock  readings  and  whose 
weight  in  EAL  was  the  sum  of  the  individual  clock  weights. 
Therefore,  the  computations  can  be  simplified  by  sumnung  the 
weights  of  the  clocks  at  each  laboratory  as  follows: 


-V Lab]^ 

WLab,=  J2^‘’ 

i=l 
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where  Nmbk  is  the  number  of  clocks  of  the  considered 
laboratory  k  and  all  the  clocks  at  a  laboratory  are  treated  as 
one  equivalent  clock  whose  reading  is  the  weighted  average 
of  individual  clock  readings.  The  formalism  of  equation  (10) 
could  still  be  applied  without  these  simplifying  assumptions,  as 
the  double  summation  would  account  for  the  100%  correlation 
of  the  time  transfer  noise  between  clocks  in  the  same  lab. 

2.2.  Extension  of  the  computation 

It  is  possible  but  not  necessary  to  apply  equation  (10)  to  every 
laboratory.  If  the  uncertainty  for  any  one  clock  is  known,  for 
example,  Xj  =  EAL  —  hj,  the  evaluation  of  uncertainty  on 
Xi  =  EAL  —  hi  may  be  obtained  by  using  the  second  equation 
in  (7)  and  applying  the  property  of  the  variance  to  obtain 

=ulj+»lj+^^xj,xu-  (12) 

The  last  term  is  the  covariance  of  the  measures  xtj  and 
the  quantity  xj  =  EAL  —  hj.  Since  all  the  clocks  are  included 
in  EAL,  the  measure  xtj  will  produce  a  nonzero  covariance 
term  by  coupling  to  the  same  measures  xtj,  which  enter  into 
the  EAL  definition  as  many  times  as  there  are  clocks  inside  or 
behind  the  laboratory  /,  obtaining 

bteqlabi 

=  -  E  (13) 

£=1 

where  Neqiabi  is  the  equivalent  number  of  clocks  in  laboratory  i , 
also  including  the  clock  external  to  that  laboratory,  but  which 
are  connected  to  UTC  through  laboratory  i.  The  equivalent 
weight  Weq^  can  be  defined  as  the  sum  of  the  Neqiabi  dock 
weights.  This  will  be  explained  further  in  the  next  section, 
where  examples  are  given.  In  this  case,  equation  (12)  becomes 

«?,  =  -  2Mx,  ,  Weq, .  (14) 

Equation  (14)  gives  the  uncertainty  of  each  clock  H/ ,  given 
the  uncertainty  of  any  clock  and  the  uncertainties  of  the 
chains  of  measures  linking  clock  H/  to  clock  Hj . 

2.3.  Correlations  and  anticorrelations 

The  BIPM  provides  the  link  uncertainties  in  Circular  T 
(section  6),  but  every  time  a  link  Xij  appears  in  a  multiple 
link  such  as  xt^k  =  +  ^j,k,  it  is  necessary  to  estimate  the 

correlations  with  respect  to  xtj,  and  this  is  the  meaning  of  the 
last  covariance  term  in  equation  (10). 

A  more  refined  evaluation  should  take  into  account  the 
different  correlation  properties  of  the  several  links,  including 
TWSTET,  GPS/Glonass  common  view  (CV)  and  Melting  Pot 
(also  termed  All-in- View)  [10,  11].  Eor  example,  in  some 
cases  the  uncertainty  is  dominated  by  calibration  uncertainties 
between  laboratory  i  and  laboratory  7,  while  in  other  cases  the 
link  uncertainty  reported  in  BIPM  Circular  T  is  mostly  due  to 
the  noise  of  the  receivers  at  both  sites. 

Correlations  will  always  occur  in  situations  where  the 
same  receiver  or  system  is  used  to  link  two  different  external 
laboratories.  The  analysis  of  these  effects  requires  more  details 
than  are  readily  available  for  correlation  of  the  links.  Eurther 


study  is  in  progress  and  a  preliminary  evaluation  indicates  that 
the  final  uncertainty  would  change  to  a  small  extent.  Some 
details  are  reported  in  the  last  section  devoted  to  future  work. 

In  this  work,  the  noise  affecting  different  nonoverlapping 
links  is  assumed  to  be  uncorrelated.  Correlations  can  and  will 
appear  when  the  same  intermediate  link  appears  in  multiple 
links,  as  will  be  shown  in  the  example  below. 

2.4.  Examples  of  application  of  the  method 

Some  examples  are  provided  in  order  to  illustrate  the  applied 
procedure.  In  each  example,  the  uncertainty  of  each  clock 
with  respect  to  UTC  increases  as  the  noise  of  the  laboratory 
k  link  increases,  but  also  decreases  as  the  weight  of  the  clock 
increases. 

Case  Al.  Let  us  consider  the  simple  case  of  only  two 
clocks  (labelled  1  and  2)  forming  UTC  and  maintained  in  the 
same  laboratory  k.  External  links  are  not  needed.  The  case  is 
depicted  in  the  sketch  below. 

k 

The  ensemble  time  UTC  would  be  computed  with  these 
two  clocks  inside  the  laboratory  k.  Suppose  clock  Hi  realizes 
UTC(k).  The  uncertainty  of  Xk  =  [UTC  —  UTC{k)f  by  means 
of  equation  (10),  would  be 

^[UTC-UTC(k)]  —  ^luTC-hi]  —  ^2^'xiX 

The  uncertainty  of  the  measure  xi  2  is  negligible  because 
the  clocks  are  within  the  same  laboratory.  Thus 

^ [UTC -UTC (k)]  ^ 

Case  A2.  Consider  two  clocks  (labelled  1  and  2),  but 
maintained  in  two  different  laboratories  (k  and  /,  respectively) 
connected  by  one  measurement  link  Xkj. 

k  I 


1  2 

In  this  case,  the  uncertainty  of  the  measurement  link  would 
be  non-negligible  and  the  final  uncertainty  for  UTC(k),  realized 
in  laboratory  k  with  clock  Hi ,  would  be 

^[UTC-UTCik)]  — 

Equation  (17)  shows  that  increasing  the  weight  of  clock 
Hi  from  laboratory  k  decreases  the  uncertainty  of  [UTC  — 
UTC{k)l 

Case  Bl.  Consider  three  clocks  (labelled  1,  2  and  3), 
maintained  in  two  different  laboratories  {k  and  /)  connected 
by  a  measurement  link  x^j  according  to  the  sketch  below. 
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k  I 


1  2  3 


For  site  k,  with  UTC(^)  realized  by  clock  Hi,  the 
uncertainty  on  [UTC  —  UTC(k)]  would  be 

^[UTC-UTC{k)]  —  '^^2^3i^ix2,k,X3,k) 

=  (W2  +  =  (1  -  (18) 

since  =  U(x2,k,x3,k)  =  because  clocks  H2  and 

H3  are  in  the  same  laboratory. 

Case  B2.  We  compute  the  uncertainty  for  UTC(/)  realized 
with  clock  H2  inside  laboratory  /  for  a  system  of  three  clocks 
(labelled  1,  2  and  3)  in  three  different  laboratories  (k,  I  and  m) 
connected  according  to  the  sketch  below. 


/ 

m 


If  the  links  between  laboratories  (/,  k)  and  laboratories 
(/,  m)  are  uncorrelated,  then  the  double  summation  in 
equation  (10)  is  zero  and 


Assuming  the  two  links  Xk,i  and  xim  are  uncorrelated,  one 
obtains  from  equation  (10) 

ulvTC-UTCik)]  =  +  «L,,]  +  2u;2W3Mu,^, 

=  (W2  +  (21) 

where  the  correlation  between  x^j  and  Xk^m  is  due  to  the 
common  path  (/,  k): 

^iXl,k,Xm,k)  —  ^Xlk’ 

This  example  shows  that,  for  uncorrelated  links,  the 
weight  of  laboratories  behind  a  pivot  laboratory  is  added  to  the 
weight  of  the  pivot  laboratory  itself  and  this  is  the  motivation 
behind  the  definition  of  ‘equivalent  weight’  Weq^  as  the  weight 
of  the  pivot  laboratory  i  plus  the  weights  of  the  clocks  behind 
that  pivot  laboratory,  which  was  introduced  in  equation  (13). 
The  last  term  in  equation  (21)  takes  into  account  the  noise  from 
the  pivot  (laboratory  1)  to  the  remote  laboratory  m,  weighted 
by  the  weight  of  clocks  in  laboratory  m. 

In  the  case  of  links  correlated  through  site-based  noise, 
with  the  same  assumptions  leading  to  equation  (20)  and 
assuming  that  the  noise  of  the  receiver  in  laboratory  I  is 
cancelled  while  connecting  laboratory  m  to  laboratory  k  with 
contemporaneous  measurements,  the  same  equations  would 
apply,  and 

^[UTC-UTC{k)]  —  (^2  +  +  ^3^^xi,k’ 

This  would  be  expected  for  links  dominated  by  site-based 
noise,  because  the  site-based  noise  of  each  site  combines  in 
such  a  way  that  all  possible  combinations  of  links  yield  the 
same  results. 

2.5.  Formulae  for  the  UTC  network 


^\uTC-UTC{l)]  —  ^i^lk,i  ^3^lm,r 

If  the  links  between  laboratories  (/,  k)  and  laboratories 
(/,  m)  are  correlated  through  site-based  noise  in  the  receiver 
in  laboratory  /,  but  the  total  noise  in  both  links  is  the  same 
and  equally  contributed  by  the  two  laboratories:  U(^xik,xi^)  = 
0.5  =  0.5  and 

U[UTC-UTC{1)]  =  “'imL  +  +  2wiW)3Mto.j,;c,,„) 

=  (w^  +  W1W3  +  (20) 

Comparing  equation  (19)  with  equation  (20),  the  effect 
of  considering  partial  correlations  between  the  link  noises 
becomes  evident. 

Case  B3.  The  same  configuration  as  Case  B2,  except  that 
the  uncertainty  for  laboratory  k  is  to  be  evaluated.  The  clocks 
in  laboratory  m  are  connected  to  lab  k  through  the  pivot  lab  / 

withX;t,m  =  Xk,l+Xi,m- 

I 

m 


Starting  from  the  examples  above,  the  uncertainty  correspond¬ 
ing  to  the  network  of  links  currently  used  for  the  computation 
of  UTC  as  reported  in  figure  1  can  be  evaluated. 

Figure  1  shows  that  USNO,  NIST  and  NICT^ ,  and  NTSC^ 
act  as  ‘intermediate’  pivots  and  the  central  role  of  the  PTB. 

Since  the  PTB  plays  a  central  role,  equation  (10)  can 
be  used  to  obtain  the  uncertainty  for  [UTC  —  UTC  (PTB)] 
considering  the  equivalent  weights  Weqi  for  the  intermediate 
pivot  laboratories  USNO,  NIST,  NICT  and  NTSC.  The  final 
expression  turns  out  to  be 


u 


[t/rc-f/rc(PTB)] 


Vi+l 

i=2 


UTC  (PTB) 


(Vi+V2)+1 

-  E 

k=Ni+2 


W^u^ 

k  UTC  (USNO) 


(Vi+V2+V3)+1 

k=(Ni+N2)+2 


k  ^Vyt,UTC(NIST) 


(Vi+V2+V3+V4)+1 

+  y]  wiu^ 

k={Ni+N2+N3)+2 


k  ^•^yt,UTC(NICT) 


+W  u  +  W  u 

^^USNO  VuTC(PTB),UTC(USNO)  ^^NIST  •^UTC(PTB),UTC(NIST) 

+  W^ 

eqmcT  •^utc(ptb),utc(nict)  c^ntsc  vutc(ptb),utc(ntsc) 

(24) 

JAl^l  VuTC(NTSC),UTC(JACT)  ^  ^ 

Here  Ai  =  25  is  the  number  oflaboratories  directly  linked 
to  PTB  with  one  single  link  (lEN,  AOS,  SP,  LT,  etc)  excluding 


^  NICT  is  the  new  name  for  CLR. 

^  All  the  acronyms  appearing  in  this  paper  are  in  agreement  with  the 
list  published  by  the  BIPM  Time  Section  at  http://www.bipm.org/jsp/en/ 
TimeFtp.jsp?TypePub=publication. 
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Figure  1.  International  time  links  used  for  computation  of  UTC  (July  2004). 


the  four  laboratories  acting  as  intermediate  pivots;  N2  =  ^ 
is  the  number  of  laboratories  linked  to  USNO  (namely  NRC, 
CNMP,  APL,  ONE  A);  A3  =  4  is  the  number  of  laboratories 
linked  to  NIST  (namely  TCC,  ONRJ,  CNM,  IGMA);  A4  =  13 
is  the  number  of  laboratories  linked  to  NICT  (SG,  TL,  NAO, 
MSL,  etc)  and  Weq^  are  equivalent  weights  of  laboratories  as 
introduced  in  equation  (13). 

The  equivalent  weight  of  the  USNO  is,  for  example, 

^e^usNO  —  ^NRC  +  WcNMP  +  E^APL  +  E^ONBA  +  E^USNO- 

Similar  expressions  hold  for  the  other  intermediate  pivots, 
NIST,  NICT  and  NTSC  and  W^ab  was  defined  in  equation  (1 1). 

Once  the  uncertainty  of  [UTC  —  UTC  (PTE)]  is  evaluated, 
the  expression  (14)  can  be  used  to  compute  the  uncertainty  of 
[UTC  —  UTCik)]  for  every  other  laboratory  k. 

3.  A  Monte  Carlo  simulation 

The  analytical  approach  presented  in  the  previous  section  was 
tested  by  a  Monte  Carlo  simulation  using  the  EIPM’s  software 
ALGOS,  on  the  full  set  of  clock  and  time  transfer  data  used  to 
generate  Circular  T published  in  July  2004.  For  that  month,  the 
pivot  laboratories  were  PTE,  USNO,  NIST,  NICT  and  NTSC. 

Since  the  source  of  uncertainty  in  [UTC  —  UTC{k)]  is 
the  links  whose  uncertainty  Uxi  j  is  listed  in  Circular  T,  a 
simulation  was  performed  assuming  that  every  link  measure  is 
described  by  a  random  variable  with  a  Gaussian  distribution. 


mean  value  equal  to  the  obtained  measurement  value  a  and 
standard  deviation  equal  to  u^^  ■ . 

The  Monte  Carlo  simulation  consists  of  picking  up 
different  values  for  the  measure  x/ ^  from  its  statistical 
distribution  and  computing  [UTC  — UTC {k)]  with  the  different 
simulated  measure  values.  This  gives  an  indication  of  the 
variability  of  the  results  [UTC  —  UTC{k)]  due  to  the  variability 
of  the  measures  which  have  been  shown  to  be  the  only 
source  of  variability  and,  hence,  of  uncertainty. 

The  simulation  proceeded  by  first  evaluating  the  effect  on 
[UTC  —  UTCik)]  of  only  one  noisy  link  (all  the  other  links  are 
considered  with  negligible  noise);  then  the  noise  was  inserted 
link  by  link.  This  ‘experimental’  variability  was  compared 
with  the  expected  theoretical  results  coming  from  the  analytical 
estimation  (10). 

3.1.  Example  of  a  single  link 

The  effect  of  the  uncertainty  of  the  link  between  the  PTE 
and  the  LT  laboratories  using  the  data  of  July  2004  was 
evaluated.  Laboratory  LT  has  one  clock  of  percentage  weight 
mu'  =  0. 1 19;  the  squared  uncertainty  of  the  link  [UTC (LT)  — 
UTC(PTB)]  is  i^luTC(LT)-uTC(PTB)]  —  27.25  ns^  and  on  MJD 
53144  [UTC(LT)  -  UTC(FTB)]  was  found  to  be  -242.8  ns. 
A  large  number.  A,  of  values  of  [UTC(LT)  -  UTC(BTB)] 
was  generated  normally  distributed  around  the  central  value 
[UTC (LT)  -  UTC(BTB)]  =  -242.8  ns  and  with  a  variance  of 
^[c/rc(LT)-c/rc(PTB)]  =  27.25  ns  . 
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[t/rc- t/rC(PTB)]/ns 


Figure  2.  Theoretical  and  experimental  distribution  of  the  values 
[UTC  —  C/rC (PTB)]  considering  only  one  link  with  noise. 


[t/rC-t/rC(LT)]/ns 


Figure  3.  Theoretical  and  experimental  distribution  of  the  values 
[UTC  —  C/rC(LT)]  considering  only  one  link  with  noise. 


Using  equation  (10),  with  only  the  link  PTB/LT  corrupted 
with  noise 

2  _  2  2 
^[f/rC-f/rC(PTB)]  —  ^LT^[C/rC(LT)-f/rC(PTB)] 

119\^ 

^ -  X  (27.25  ns^)  =  3.85  x  10“^  nsl 

100  / 


Thus 

W[[/rc-t/rc(PTB)]  =  0.0062  ns, 


and  from  equation  (14) 

W[[/rc-t/rc(LT)]  =  5.21  ns. 


Alternatively,  using  equation  (10)  directly  referenced  to 
laboratory  LT,  which  corresponds  to  having  the  total  clock 
weight  in  PTB  with  the  exception  of  the  LT  clock,  the  same 
theoretical  result  was  obtained: 

/  0.119\  ! - T 

W[f/rc-f/rc(LT)]  ~  I  ^ - I  ^  V27.25ns^  =  5.21ns. 

Using  N  =  10000  simulated  measurement  values,  a 
normal  distribution  of  values  [UTC— UTC  (PTB)]  was  obtained 
with  a  mean  value  equal  to  —  1 6.7087  ns  and  standard  deviation 
equal  to  0.0063  ns.  The  value  [UTC  -  UrC(PTB)]  for  that 
date  published  in  Circular  T  was  equal  to  —16.7  ns  and 
the  standard  uncertainty  from  the  above  computation  was 
i^[t/rc-LTC(PTB)]  =  0.0062  ns.  Therefore,  the  mean  value 
and  the  standard  deviation  of  [UTC  —  UTC(PTB)]  obtained 
by  simulating  the  noisy  link  PTB/LT  corresponded  to  the 
published  value  and  to  the  expected  theoretical  uncertainty. 

For  the  laboratory  LT  a  normal  distribution  of  mean  value 
equal  to  226.161  ns  was  obtained,  with  a  standard  deviation 
of  5.2ns.  The  published  value  in  Circular  T  was  [UTC  — 
UrC(LT)]  =  226.1  ns  and  the  standard  uncertainty  from  the 
above  theoretical  computation  was  u\jjtc-utc{\.t)^  =  5.21  ns. 

The  results  are  depicted  in  figures  2-4,  where  the 
histogram  reports  the  results  obtained  from  the  Monte  Carlo 
simulation,  while  the  solid  line  (in  red)  represents  the  expected 
statistical  distribution  of  the  values  according  to  the  theoretical 
analysis. 

As  can  be  seen,  the  analytical  and  simulation  results  agree. 


Figure  4.  Theoretical  and  experimental  distribution  of  the  values  of 
[UTC  -  UTC(mST)l 


3.2.  Complete  system 

The  Monte  Carlo  results  are  reported  here  for  the  complete 
system  of  links,  using  A  =  20  000  simulated  data  for  any 
link.  Figure  4  shows  the  close  agreement,  for  example, 
between  the  theoretical  and  experimental  results  for  the  values 
of  [UTC  —  UrC(NIST)].  In  table  1  the  analytical  results  are 
compared  to  the  simulations  assuming  that  every  link  had  the 
noise  corresponding  to  the  uncertainty  listed  in  Circular  T.  The 
agreement  between  the  two  estimation  methods  is  evident  and 
gives  confidence  to  the  analysis. 


4.  Application  to  the  computation  of  Circular  T 

The  formulae  developed  in  section  2.5  for  the  evaluation  of 
[UTC  —  UTC(k)]  uncertainties  have  been  applied  using  the 
link  uncertainties  as  listed  in  section  6  of  Circular  T.  They 
have  been  used  operationally  since  January  2005,  reporting 
the  uncertainties  in  Circular  T.  Results  for  Circular  T  no  205 
of  January  2005  are  reported  in  figure  5. 
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Table  1.  Analytical  and  numerically  estimated  uncertainties  of  all  the  UTC  participating  laboratories  considering  every  link  with  noise. 


k 

U[UTC-UTC(k)]^^^ 

k 

U[UTC- 

c/rc(A:)]/ns 

Analytical  method  Numerical  method 

Analytical  method 

Numerical  method 

AOS 

5.6 

5.6 

NIS 

20.3 

20.6 

APL 

5.7 

5.7 

NIST 

5.0 

5.0 

AUS 

7.3 

13 

NMC 

20.7 

20.8 

BEV 

5.5 

5.5 

NMIJ 

7.0 

7.0 

BIRM 

20.6 

20.6 

NPL 

5.3 

5.2 

CAO 

21.2 

21.0 

NPLI 

20.3 

20.4 

CH 

5.4 

5.3 

NRG 

15.2 

15.1 

CNM 

21.0 

21.1 

NTSC 

7.1 

7.0 

CNMP 

8.4 

8.4 

OMH 

20.2 

20.5 

CSIR 

20.3 

20.6 

ONBA 

8.9 

8.9 

DLR 

5.4 

5.4 

ONRJ 

21.2 

21.3 

DTAG 

10.6 

10.5 

OP 

5.4 

5.4 

HKO 

7.3 

13 

ORB 

5.3 

5.4 

lEN 

2.4 

2.4 

PL 

5.4 

5.3 

lEAG 

5.3 

5.3 

PTB 

1.9 

1.9 

IGMA 

20.8 

20.9 

ROA 

5.4 

5.4 

INPL 

10.9 

11.0 

SCL 

11.6 

11.6 

JATC 

21.2 

21.3 

SG 

20.5 

20.6 

JV 

20.7 

20.6 

SMU 

20.7 

20.7 

KRIS 

7.1 

12 

SP 

10.4 

10.5 

EDS 

20.3 

20.2 

SU 

6.1 

6.1 

LT 

5.6 

5.6 

TCC 

21.2 

20.9 

MSL 

20.7 

20.8 

TL 

6.7 

6.7 

NAO 

20.6 

20.7 

TP 

5.8 

5.8 

NICT 

4.4 

4.4 

UME 

25.1 

25.1 

NIM 

20.4 

20.3 

USNO 

2.3 

2.3 

NIMT 

20.7 

20.9 

VSL 

5.3 

5.3 
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Date 

2004/05  Oh  UTC 

DEC  30 

JAN  4 

JAN  9 

JAN  14 

JAN  19 

JAN  24 

JAN  2  9 

Uncertainty/ns 

Notes 

MJD 

53369 

53374 

53379 

53384 

53389 

53394 

53399 

Laboratory  k 

{UTC-UTC{k) ] /ns 

Ua 

Ub 

u 

AOS 

(Borowiec) 

-21.2 

-14.8 

-18.7 

-13.6 

-11.2 

-12.1 

-17.9 

1.6 

5.4 

5.6 

BEV 

(Wien) 

84.9 

87.5 

90.2 

90.2 

94.8 

96.9 

99.9 

1.5 

5.3 

5.5 

CH 

(Bern) 

-27.4 

-25.2 

-27.4 

-31.1 

-28.6 

-23.8 

-21.1 

0.8 

5.3 

5.4 

DLR 

(Oberpfaf f enhof en) 

1.4 

-3.1 

-7.8 

-14.1 

-19.1 

-22.3 

-31.3 

0.8 

5.4 

5.5 

lEN 

(Torino) 

-58.2 

-71.8 

-77.9 

-79.8 

-78.5 

-84.6 

-87.5 

0.7 

2.2 

2.3 

LT 

(Vilnius) 

486.0 

519.1 

547.4 

549.6 

528.4 

533.6 

542.7 

1.6 

5.3 

5.5 

NICT 

(Tokyo) 

4.5 

-0.7 

6.8 

6.1 

-1.7 

1.1 

3.9 

1.2 

4.3 

4.5 

NIST 

(Boulder) 

3.5 

2.6 

1.3 

1.5 

0.8 

-0.6 

-1.7 

0.6 

5.1 

5.1 

NPL 

(Teddington) 

8.0 

11.2 

14.6 

18.5 

22.7 

25.7 

33.7 

0.7 

2.2 

2.3 

NTSC 

(Lintong) 

17.7 

20.4 

26.7 

31.1 

27.4 

19.6 

19.1 

2.7 

6.5 

7.0 

ONBA 

(Buenos  Aires) 

-923.0 

-967.0 

-1200.8 

-1407.0 

-1666.2 

-3418.7 

-3392.1 

5.0 

7.3 

8.8 

ONRJ 

(Rio  de  Janeiro) 

155.1 

153.1 

162.1 

159.6 

161.7 

158.4 

176.0 

5.0 

20.6 

21.2 

OP 

(Paris) 

41.7 

45.8 

49.4 

43.6 

40.3 

39.8 

36.4 

0.6 

2.1 

2.2 

ORB 

(Bruxelles) 

-42.9 

-41.2 

-39.3 

-37.5 

-35.2 

-34.5 

-32.3 

0.8 

5.3 

5.4 

PL 

(Warszawa) 

-3.9 

-1.8 

-5.7 

-6.5 

-4.5 

-4.1 

-7.8 

1.5 

5.2 

5.4 

PTB 

(Braunschweig) 

-0.6 

1.2 

2.5 

3.0 

5.5 

5.9 

4.7 

0.4 

1.9 

1.9 

ROA 

(San  Fernando) 

-63.1 

-62.0 

-68.3 

-67.9 

-67.1 

-65.8 

-74.1 

0.8 

5.3 

5.4 

UME 

(Gebze-Kocaeli) 

599.8 

595.1 

593.2 

603.0 

608.3 

611.2 

612.5 

15.0 

20.1 

25.1 

USNO 

(Washington  DC) 

-4.9 

-4.7 

-5.1 

-4.6 

-3.5 

-3.5 

-2.5 

0.5 

2.2 

2.3 

VSL 

(Delft) 

-76.9 

-77.0 

-77.9 

-77.4 

-56.4 

-41.8 

-21.2 

0.7 

2.2 

2.3 

Figure  5.  A  sample  of  the  first  section  of  Circular  T  no.  205  showing  uncertainties  in  [UTC  —  UTC{k)\ 
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In  Circular  T,  the  last  column  (see  figure  5),  denoted 
by  u,  corresponds  to  the  combined  standard  uncertainty  of 
[UTC  —  UTC(k)]  for  any  laboratory  k.  The  values  ua  and 
ub  are  the  components  of  the  combined  standard  uncertainty 
u  originating  from  type  A  and  type  B  evaluations  listed  in 
section  6,  respectively.  The  following  relationship  holds: 


u  = 

It  is  noted  that  the  application  of  equations  (24)  and  (14) 
separately  to  the  type  A  and  type  B  evaluations  is  permitted 
only  when  all  the  other  quantities  involved  are  constant,  in 
particular  the  number  of  clocks  and  their  weights  are  constant; 
this  is  currently  true  only  within  each  separate  month  of  UTC 
computation. 

The  application  of  the  above  theory  to  the  actual  case  of 
[UTC  —  UTC(k)]  suggests  an  additional  experimental  test  to 
check  the  validity  of  the  obtained  estimates.  In  fact,  under  the 
assumptions  that  the  values  of  [UTC  —  UTC(k)]  are  almost 
constant  or  very  slowly  varying  in  time,  which  is  a  reasonable 
assumption  in  case  of  the  UTC(^)  time  scale  realized  with 
hydrogen  masers,  one  can  ascribe  the  observed  instability  of 
the  values  [f/rC  —  UTC (k)],  at  a  five-day  observation  interval, 
mostly  to  the  combined  effect  of  time  transfer  noise  rather  than 
to  clock  noise.  By  estimating  the  Allan  deviation  or  the  Time 
Deviation  (TDEV)  of  the  quantity  [UTC  —  UTC(k)]  one  has 
thus  a  rough  idea  of  the  combined  time  transfer  noise  of  all  the 
links  in  TAI  as  observed  from  laboratory  k  and  therefore  gets  an 
independent  estimate  of  the  ua  uncertainty  for  that  laboratory. 
A  simple  test  evaluating  the  TDEV  for  a  few  time  laboratories 
gave  reasonably  good  agreement  with  the  calculated  values  of 
ua  reported  here. 

5.  Future  extension 

To  refine  the  evaluation  of  uncertainties,  it  is  necessary  to  have 
more  details  than  are  readily  available  about  the  correlation  of 
the  links.  This  is  largely  because,  as  noted  in  the  examples, 
the  correlated  noise  between  two  sites  should  include  the 
contribution  from  site-based  uncertainties  at  sites  along  the 
path  of  links  and  because  the  correlation  between  two  links 
with  a  common  site  that  uses  the  same  equipment  for  both 
links  should  include  the  site-based  noise  contribution  from  that 
site.  This  is  difficult  to  estimate  when  calibrations  are  done  by 
links  rather  than  by  sites.  In  computations  assuming  site-based 
noise,  which  supplemented  the  link  uncertainties  reported  in 
Circular  T,  it  was  found  that  on  the  whole  the  effects  of  these 
simplifications  are  small. 

Eurther  extension  of  this  evaluation  would  be  based  on  the 
following  consideration.  If  all  time  transfers  were  achieved 
using  a  single  system  per  site  and  if  all  sources  of  noise  were 
site-based,  such  as  mostly  happens  for  Melting  Pot  GPS,  also 
termed  All-in- View  (AV)  [10, 1 1],  then  all  possible  links  would 
obey  the  following  closure  relation: 

•T/,j(0  +Xj^^(t)  +  (0  =  0.  (25) 

In  this  situation,  the  noise  of  each  site’s  time  transfer 
system  would  be  indistinguishable  from  the  noise  of  its  clocks 
and  the  dominant  uncertainty  would  be  given  by  the  site  noise. 


and  all  the  external  clocks  would  be  seen  through  that  dominant 
noise  independently  of  their  location.  If  different  satellite 
schedules  are  used  in  the  relevant  links,  uncertainties  in  time 
transfer  using  GPS  CV  are  largely,  though  not  entirely,  site- 
based.  Even  for  CV  observations  made  using  every  available 
satellite,  closure  violations  will  only  arise  if  simultaneous 
satellite  observations  are  recorded  at  only  two  of  the  three 
sites.  In  such  cases,  orbit  misestimation  and  receiver  noise 
will  contribute  uncertainties,  and  any  azimuth  or  elevation- 
dependent  asymmetries  in  the  multipath  environment  would 
cause  both  uncertainties  and  biases.  Since  calibration  is 
achieved  by  an  all-sky  sampling  that  is  systematically  different 
from  the  sky- sampling  of  CV,  a  different  multipath  will 
also  lead  to  uncertainties.  Despite  these  noise  sources,  the 
closure  relation  largely  holds  for  CV,  and  the  largest  source  of 
uncertainty  is  typically  due  to  variations  of  the  receiver  system 
that  are  common  to  all  data. 

Eor  Two-Way  Satellite  Time  and  Erequency  Transfer 
(TWSTET),  the  noise  is  again  largely  site-dependent.  Some 
closure  violations  can  occur  because  the  observations  between 
pairs  of  sites  are  typically  made  with  different  codes  and 
slightly  different  frequencies.  The  largest  source  of  closure 
errors  is  probably  due  to  the  fact  that  the  received  signals 
are  shaped  by  the  product  of  the  transmitting  and  receiving 
bandpasses,  while  the  delay  and  certain  noise  components 
such  as  the  cable-dependent  multipath  can  systematically  vary 
over  the  bandpass  [12].  While  TWSTET  closure  violations 
are  seen  at  the  1  ns  level  in  the  data  sent  to  the  BIPM,  they 
could  be  reduced  through  baseline-dependent  calibrations.  In 
the  cases  where  a  TWSTET  system  is  calibrated  with  GPS, 
the  uncertainties  in  the  calibration  are  determined  by  the 
uncertainties  in  the  GPS  calibration. 

Special  situations  arise  when  one  site  is  a  pivot  site, 
connected  to  some  sites  by  one  technique  and  other  sites  by 
a  different  technique.  By  means  of  an  illustration,  we  can 
consider  a  very  simplified  situation  in  which  every  site  is 
directly  linked  to  one  central  pivot  site,  either  by  AV  GPS 
or  TWSTET;  this  can  easily  be  generalized  to  describe  more 
complex  topologies.  We  will  assume  that  variations  between 
any  two  links  are  completely  uncorrelated  with  respect  to 
variations  at  the  link  extrenfities,  but  that  the  links  are  100% 
correlated  with  respect  to  variations  of  the  equipment  at  the 
central  pivot  site,  provided  both  links  are  by  either  TWSTET  or 
GPS.  Let  us  also  assume  a  bias  B  exists  in  the  GPS  equipment 
at  the  pivot  site.  In  this  case,  it  is  easy  to  show  that  the  bias 
would  affect  those  laboratories  k  linked  by  GPS  to  the  pivot  as 
follows: 

^[UTC-UTC{k)]  =  (1  —  Wg)5 

where  Wq  is  the  sum  of  all  the  weights  of  the  laboratories 
linked  to  the  pivot  site  by  GPS. 

The  bias  would  affect  the  pivot  laboratory  and  those  linked 
to  it  by  TWSTET  as  follows: 

^[UTC-UTC(k)]  =  —WqB. 

Under  the  normal  circumstances  described  in  this  paper, 
the  existence  of  any  biases  would  not  carry  any  significant 
statistical  implications,  as  they  would  be  directly  related  to 
the  tabulated  uncertainties  in  the  links  themselves.  However, 
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the  above  equations  illustrate  the  dependence  of  TAI  and  UTC 
upon  the  equipment  at  any  central  pivot,  which  may  not  always 
follow  the  Gaussian  behaviour  assumed  here,  particularly  in 
the  case  of  equipment  failure. 

6.  Conclusions 

This  paper  presents  a  study  of  the  determination  of  the 
uncertainties  in  {UTC  —  UTC{k)].  An  analytical  solution 
is  derived  from  the  law  of  the  propagation  of  uncertainty, 
taking  into  account  that  leap  seconds  and  deterministic 
frequency  steering  of  EAL  do  not  affect  these  uncertainties. 
The  analytical  results  were  verified  through  Monte  Carlo 
simulations  using  the  software  used  to  generate  UTC,  and 
agreement  was  found,  giving  confidence  in  the  analytical 
estimation. 

A  more  detailed  analysis  is  in  progress,  including  full 
inclusion  of  site-based  noise  contributions,  all  available 
calibration  information,  more  details  for  the  correlation 
between  the  links,  methods  for  optimizing  the  link  structure, 
given  uncertainty  information,  non-Gaussian  behaviour  and 
different  correlation  properties  of  uncertainties  due  to 
calibration  or  due  to  random  noise. 
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